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Work  Statement 


1.  Determine  experimentally  the  conditions  under  which  free-floating  nitrogen  and 
helium  plasmas  can  be  initiated  and  sustained  in  the  TEq^  mode  microwave  resonant 
cavity.  Measure  maximum  obtainable  gas  pressure  and  temperature  and  coupling 
efficiency  as  a  function  of  input  power  and  gas  flow  rate.  Spectroscopic 
measurements  of  temperature  and  e;ectron  density  will  not  assume  complete  LTE. 
Compare  results  to  those  obtained  for  the  TMoi  mode. 

2.  Examine  the  effect  of  various  plasma  stabilization  schemes  in  a  microwave  resonant 
cavity  on  maximum  obtainable  gas  pressure  and  temperature  and  coupling  efficiency. 
Quantify  the  losses  caused  by  the  introduction  of  the  plasma  stabilization  schemes. 

3.  Modify  the  numerical  model  of  the  microwave-heated  propagating  plasma  to 
incorporate  the  effects  of  energy  transfer  to  the  cold  upstream  gas  by  resonant 
radiation.  Compare  numerical  results  to  experimentally  measured  values. 

4.  Examine  the  effects  of  various  techniques  to  spatially  stabilize  experimentally 
generated  propagating  plasmas.  Measure  plasma  and  gas  temperature.and  coupling 
efficiencies  as  a  function  of  input  power  for  the  stabilized  plasmas. 
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1-P  vStatus  of  Research  Effort 

Experiments  in  the  previous  grant  year  with  free-floating  resonant  cavity  plasmas 
indicated  that  the  maximum  input  power  was  restricted  to  below  approximately  500  W, 
above  which  the  plasma  was  observed  to  move  towards  the  wall  of  the  quartz  vessel, 
resulting  in  rapid  heating  of  the  quartz.  It  is  believed  that  this  instability  is  due  to  the 
asymmetric  introduction  of  microwave  power  to  the  cavity,  which  results  in  a  perturbation 
of  the  electric  field  pattern  from  the  axisymmetric  pattern  which  would  otherwise  be 
expected.  In  order  to  overcome  this  instability  and  allow  the  maximum  available  power  to 
be  input  to  the  resonant  cavity  (2.5;kW)  two  fluid  dynamic  stabilization  schemes,  bluff  body 
and  swirling  flow,  were  examined  this  year.  Both  these  schemes  have  been  used  extensively 
in  combustion  systems  and  to  some  extent  in  previous  electrothermal  propulsion  systems 
involving  plasmas.  The  bluff  body  creates  a  recirculation  region  with  a  decreased  pressure 
and  increased  residence  time  while  enhancing  mixing  in  the  region  behind  the  bluff  body 
while  swirling  flows  also  create  a  region  of  decreased  pressure  and  a  recirculating  flow 
without  the  fluid  dynamic  drag  associated  with  the  bluff  body  and  the  heat  transfer  to  the 
bluff  body. 

The  plasma  containment  vessel  used  in  these  experiments  was  a  quartz  hemisphere 
of  75  mm  diameter  with  a  boron  nitride  flange  at  the  downstream  end  containing  a  choked 
orifice.  The  bluff  body  and  flow  s^virler  were  also  fabricated  from  boron  nitride.  Both 
helium  and  nitrogen  plasmas  were  investigated  and  the  resonant  cavity  was  operated  in. the 
TMqii  mode  (compared  to  the  TMqis  mode  used  for  the  free-floating  plasmas),  which  gave 
a  region  of  maximum  electric  field  strength  in  the  region  between.the  flow  stabilizer  and  the 
choked  orifice.  Maximum  gas  pressure. available  with  the  experimental  system  was  500  kPa. 
(5  atmospheres). 

It  was  found  that  bluff-body  stabilization  proved  more  effective  in  allowing  operation 
at  higher  input  microwave  power  levels  and  gas  pressures  than  the  flow  swirler,  although  the 
poor  performance  of  the  flow  swirler  may  be  due  to  its  installation  upstream  of  the 
hemispherical  quartz  vessel  and  the  resulting  expansion  of  the  flow  as  it  enters  the 
hemisphere.  Thus  the  majority  of  the  experiments  were  carried  out  using  the  bluff  body. 
It  was  also  found  that  with  the  bluff  body  nitrogen  plasmas  could  only  be  sustained  at 
powers  less  than  700  W  and  gas  pressures  less  than  170  kPa  before  becoming  unstable  and 
moving  to  the  wall  of  the  quartz  vessel,  whereas  helium  plasmas  could  be  sustained  up  to 
the  maximum  available  microwave  power  of  2;5  kW  and  gas  pressure  of  500  kPa.  The 
majority  of  the  measurements  were  thus  made  with  helium  plasmas. 

Fi^re  1  shows  the.  ability  of  the  bluff  body  to  raise  the  maximum  microwave  power 
and  gas  pressure  while  keeping  a  helium  plasma  centered  in  the  resonant  cavity.  The  points 
plotted  for  the  TMqjj  mode  are  for  the  free-floating  plasma  where  higher  input  powers 
resulted  in  the  plasma  moving  to  the  wall  of  the  resonant  cavity.  The  points  plotted  for  the 
TMoii  mode  are  for  the  bluff-body  stabilized  plasma.  It  can  been  that  for  low  power  levels 
the  results  for  the  two  modes  coincide  while  the  bluff  body  allows  operation  at  much  higher 
power  levels  and  gas  pressures.  The  maximum  gas  pressures  of  500  kPa  at  1000  and  1500 
W  do  not  indicate  that  the  plasma  cannot  be  sustained  at  higher  pressures  but  is  only  the 
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limitation  of  the  experimental  systein.  At  these  points  no  increase  of  reflected  power  from 
the  cavity  was  observed,  an  increase  which  serves  as  an  indication  of  approaching  the 
extinction  limit. 

The  coupling  efficiency  is  a  measure  of  the  percentags  of  the  input  microwave  power 
which  is  absorbed  by  the  cavity  and  is  the  ratio  of  the  input  power  minus  the  reflected 
power  divided  by  the  input  power.  Figure  2  plots  the  coupling  efficiency  as  a  function  of 
gas  pressure  for  bluff-body  stabilized  helium  plasmas  for  input  powers  from  1.25  to  2  kW. 
It  can  be  seen  that  the  coupling  efficiency  increases  with  pressure,  with  higher  input  powers 
requiring  higher  pressures  to  obtain  their  maximum  coupling  efficiency.  In  those  cases 
where  the  experimental  system  permitted  a  sufficiently  high  pressure  to  reach  the  maximum, 
coupling  efficiencies  close  to  100%  where  obtained.  Even  with  the  low  powers  and 
pressures  which  allowed  a  stable  nitrogen  plasma,  coupling  efficiencies  approaching  100% 
were  obtained.  These  results  are  much  higher  than  the  80%  coupling  efficiencies  obtained 
with  the  free-floating  helium  and  nitrogen  plasmas  operating  in  the  TM012  mode  because  the 
bluff-body  stabilization  permitted  more  optimal  tuning  of  the  cavity  \vithout  the  plasma 
moving  to  the  wall  of  the  quartz  vessel,  thus  indicating  another  advantage  of  the  fluid 
dynamic  stabilization  provided  by  the  bluff  body. 

The  overall  efficient  of  the  microwave  electrothermal  thruster  is  the  ratio  of  the 
power  of  the  gas  exiting  through  the  choked  orifice  to  the  input  microwave  power.  The 
power  of  the  gas  exiting  through  the  choked  orifice  was  determined  by  the  mean  stagnation 
temperature  of  the  gas  and  the  mass  flow  rate.  The  mean  stagnation  temperature  of  the  gas 
was  determined  by  measuring  the  stagnation  pressure  and  rnass  flow  rate  through  the  quartz 
vessel  and  knowing  the  area  of  the  choked  orifice.  Figure  3  plots  the  overall  efficiency  of 
bluff-body  stabilized  helium  plasmas  as  a  function  of  specific  power  (microwave  power 
absorbed  by  the  plasma  divided  by  the  mass  flow)  for  input  powers  from  500  W  to  2  kW. 
The  processes  which  cause  an  overall  efficiencty  of  less  than  100%  are  conductive  losses  to 
the  bluff  body,  orifice  plate  and  quartz  walls,  reflected  power  from  the  cavity,  ohmic  losses 
in  the  cavity  walls  and  radiative  losses  from  the  plasma,  which  will  be  absorbed  by  either 
the  quartz  vessel  or  the  surrounding  cavity.  For  a  given  input  power  the  specific  power  was 
increased  by  decreasing  the  chamber  pressure  and  therefore  the  mass  flow  rate  through  the 
choked  orifice.  It  can  be  seen  that  the  overall  efficiency  decreases  with  decreasing  pressure 
(and  therefore  mass  flow  rate)  and  increasing  power.  This  may  be  explained  by  considering 
the  physical  behavior  of  the  plasma  as  these  parameters  are  varied.  Increasing  the  pressure 
reduced  the  plasma  diameter.  In  the  region  where  the  plasma  is  stably  located  downstream 
of  the  bluff  body,  this  reduction  in  plasma  diameter  has  the  effect  of  reducing  the  heat  load 
to  the  quartz  vessel  and  thus  increasing  the  overall  efficiency.  Increasing  the  input  power 
has  the  effect  of  increasing  the  plasma  diameter  (as  well  as  the  plasma  length),  increasing 
the  conductive  losses  to  the  bluff  body,  orifice  plate  and  quartz  vessel,  therefore  reducing 
the  overall  efficiency. 

Assuming  an  isentropic  expansion  to  vacuum  conditions,  the  specific  impulse  of  the 
rhicrowave  electrothermal  thmster  was  calculated.  It  was  found  that  the  specific  impulse 
increased  linearly  with  specific  power,  reaching  a  maximum  value  of  543  s  for  the  helium 
plasma.  The  use  of  hydrogen  propellant  under  the  same  operating  conditions  would 


produce  a  specific  impulse  of  1086  s  due  to  the;  decrease  of  molecular  weight  (assuming 
complete  dissociation,  which  Would  be  expected  at  the  plasma  temperatures  achieved)  by 
a  factor  of  four. 

Spectroscopic  measurements  were  made  of  the  helium  pl^ma  temperature  as  a 
function  of  the  operating  conditions.  The  spectroscopic  system  consisted  of  a  Czemy-Tumer 
half  meter  spectrometer  and  an  electronically  tunable  Fabry-Perot  etalon.  Electron 
temperatures  were  measured  using  the  absolute  continuum  method.  Calibration  of  the 
system  was  accomplished  by  using  a  calibrated  tungsten  strip  lamp  and  the  Saha  equation 
was  used  to  obtain  the  electron  number  density  as  a  function  of  temperature  (the  electron 
number  density  being  a  required  input  parameter  to  determine  the  electron  temperature  via 
this  method).  Measurements  of  the  Doppler  broadened  line  width  were  used  to  determine 
the  heavy  particle  (ions  and  neutrals)  temperature.  Radial  profiles  of  electron  temperature 
were  measured  using  an  Abel  inversion.  It  was  found  that  the  electron  temperature  at  the 
center  of  the  plasma  was  insensitive  to  variations  in  input  power  and  gas  pressure,  varying 
from  11,880  to  12,170  K  as  the  input  power  was  varied  from  500  to  1000  W  and  the  gas 
pressure  from  120  to  300  kPa  (absolute).  This  insensitivity  corresponds  to  results  published 
in  previous  literature  for  microwave-heated  argon  and  hydrogen  plasmas.  Radial  profiles 
of  the  electron  temperature  from  the  center  to  the  edge  of  the  plasma  were  found  to  be  flat, 
decreasing  approximately  10%  from  the  center  to  the  edge.  Measurements  of  the  heavy 
pa,rticle  temperature  were  more  sensitive  to  the  light  intensity  measurements  and  therefore 
produced  greater  scatter  in  the  results,  however  the  data  scatter  bounded  the  measured 
electron  temperatures  with  the  conclusion  that  the  electron  and  heavy  particle  temperatures 
are  approximately  equal. 

Because  the  resonant  cavity  and  waveguide  experiments  could  not  be  conducted 
simultaneously  since  they  both  required  the  same  microwave  power  source,  the  final 
experiments  with  the  propagating  microwave  waveguide  plasmas  could  not  be  completed  by 
the  end  of  the  grant  period.  However  these  experiments  are  continuing  under  support  by 
the  NASA  Propulsion  Engineering  Research  Center  at  Penn  State.  Previous  experiements 
with  bluff-body  stabilized  heliuni  and  nitrogen  waveguide  plasmas  used  a  rectangular 
waveguide  which  produced  a  plasma  which  extended  from  one  wall  of  the  waveguide  to  the 
other  due  to  the  electric  field  pattern  propagating  through  the  waveguide.  This  plasma 
geometry  is  not  desirable  from  a  propulsion  standpoint  because  of  the  contact  between  the 
plasma  and  the  walls  of  the  wave^ide  and  the  associated  erosion  and  heat  transfer  losses. 
A  cylindrical  waveguide  was  designed  and  manufactured  which  propagates  an  electric  field 
where  the  field  maxima  are  in  the  interior  of  the  waveguide,  producing  a  plasma  which  is 
centrally  located,  not  making  contact  with  the  walls  of  the  cavity.  A  bluff  body  is  once.again 
used  to  provide  fluid  dynamic  stabilization  of  the  plasma  with  the  flow  exiting  through  a 
choked  orifice.  Experiments  are  currently  underway  with  helium  plasmas  to  measure 
stability  boundaries,  coupling  efficiencies,  overall  efficiencies  and  electron  and  heavy  particle 
temperatures  as  a  function  of  input  microwave  power  and  gas  pressure. 
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Efficiency  (%) 


Fig.  2  Coupling  efficiencies  for  bluff-body  stabilized  helium  plasm^  as  a  function  of 
pressure  for  input  powers^  from  125  to  2  kW.  Higher  powers  require:higher 
pressures  for  optim^  coupling  efficiency.. 
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Efficiency  (7.) 


Specific  Power  (MJ/kg) 


Fig.  3  Calculated  overall  efficiencies  versus  specific-power  for  bluffcbody  stabilized 

heliuni  pl^m^  shdvdngia  decre^c/in  overall  effidency  down  tn  a^mim'tnnm 
of  45%  with  an  increase  in  specific  power. 
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2.0  DESCRIPnOJ  OF  THE  (XMTOTATIONftL  MODEL 

TSie  ccnpitational  model  for  micaxwave  thermal  propulsion  is  based 
v?xai  the  simultaneous  numerical  solution  of  the  Maxwell  and  Navier- 
Stokes  equations.  These  equations  are  strongly  cot^jled  to  each  other 
throu^  the  electrical  ccxiductivity  of  the  gas  vMch  acts  as  a  source 
term  in  the  energy  equatiai  to  absorb  the  micrawaves,  and  also  as  a 
parameter  in  the  Maxwell  equaticxis  that  determine  the  electronagnetic 
field  distribution.  Althou^  previous  analyses  of  radio  frequency 
plasmas  [1-3]  have  solved  the  Maxwell  equations  in  the  frequenqi'^ 
donain,  the  present  analysis  solves  the  conplete  time-di^)end^t 
equations  in  the  physical  donain.  ■  The  harmonic  analysis  is  suitable 
vhen  the  fundamental  mode  is  of  interest,  but  its  convergence  is 
seriously  impaired  for  hi^er  modes  such  as  the  TM012  and  the  TMOll 
modes  that  are  of  interest  here.  In  particular,  an  explicit  time- 
dqpendent  ^proach  similar  to  Yee's  [4]  is  chosen.  The  EMS  values  of 
•Uie  fields  required  in  the  Navier-Stdkes  equations  are  determined  by 
averaging  the  instantaneous  E-M  field  over  the  wave  period.  To  find 
this  value,  the  E-M  equations  are  solved  over  one  (or  several)  wave 
period(s)  at  each  time  step  of  the  wavier-Stokes  equations.  Because 
the  explicit  scheme  i:ised  for  the  Maxv^l  equations  is  extremely  rapid, 
this  procedure  causes  less  than  a  factor  of  two  increase  in  the  CEU 
time  required  as  oonpared  to  the  Navier-Stokes  equations  alone. 

Parametric  soluticMis  have  been  obtained  for  three  basic 
geometrical  configurations  for  micrcwave  absorption,  the  free-floating 
plasma  in  a  resonant  cavity,  the  bluff-body  stabilized  plasma  in  a 
resonant  cavity,  and  the  bluff-bocfy  stabilized  plasma  in  a  waveguide. 
These  three  geometrical  configurations  are  identical  to  the  three 
ejqjerimental  configurations  that  have  been  tested,  and  comparisons 
with  ei^jeriment  are  given  for  all  three  geometries.  The  same  analysis 
is  used  for  all  three  configurations.  This  model  and  the  numerical 
solvition  technique  is  outlined  in  this  sectic^,  follcwed  by  a  summary 
of  the  results  in  the  following  section. 

2.1  Geometric  Oonficrurations 

The  three  geometric  configurations  considered  are  given  in 
Figs.  1  and  2.  Figure  la  shews  the  geometry  for  the  free-floating 
resonant  cavity  plasma,  \^iile  Fig.  lb  r^resents  the  bluff-botfy 
stabilized  resonant  cavity  plasma  and  Fig.  2  d^icts  the  bluff-bocty 
stabilized,  waveguide  heated  plasma.  The  free-floating  resonant 
cavity  plasma  in  Fig.  la  consists  of  a  sphere-cylinder  tube 
arrangement  that  passes  threu^  the  micrt»ave  cavity.  The  tube 
diameter  is  2.5  on  while  the  ^here  diameter  is  10  cm.  In  general,  a 
TM012  node  is  set  rp  in  the  cavity,  and  is  used  to  sustain  a  micrcwave 
plasma  in  the  spherical  section  of  the  cavity.  To  understand  the 
details  of  the  fluid  c^namic  interactions  in  this  geemetry  better,  we 
have  cdso  cemputed  solutions  for  a  ^here-cylinder  combination  with  a 
tube  diameter  of  7.5  cm,  and  a  strai<^t  tube  configuration  in  vMch 
the  entire  flcwfield  is  cylindrical  in  nature  (no  broadened  spherical 
region) ,  with  a  diameter  of  10  cm.  Examples  of  these  special 
geometries  are  given  later.  In  all  cases,  the  flew  is  from  top  to 
bottom,  with  gravitational  forces  acting  in  the  cJcvaiward  direction. 
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The  gecraetrical  details  of  the  flew  tube  and  the  itdcxcwave  cavity 
for  the  bluff-4xxty  stabilized  rescanant  cavity  plasma  are  described  in 
Fig.  lb.  Hie  gas  agpin  flows  from  top  to  bottom  in  a  cylindrical  tube 
with  a  diameter  of  4  cm,  with  a  wider  bulb  at  the  bottom.  The  bulb 
ejqpands  out  to  a  diameter  of  8  cm,  and  the  gas  exits  from  this  bulb 
throu^  a  small,  choked  hole  in  the  bottem.  Hie  bluff  bo(^  is  a  small 
triangular^-shc^ied  (non-ccmducting)  boefy  with  a  base  diameter  of  2  an. 
In  the  experiments,  the  bluff  bocty  can  be  moved  axially  and  dewn  to 
find  the  most  "effective”  locatiai  for  the  plasma.  Hie  standing 
electromagnetic  field  pattern  in  this  ocnfiguratiai  is  a  HMOll  mode. 
Hie  plasma  again  forms  at  one  of  the  peaks  in  the  standing  mode.  A 
second  alternate  geemetry  is  also  used  for  broadened  understanding  in 
this  geometrical  ccxifiguration. 

Hie  waveguide  configuration  is  a  strai^t  tube  configuration 
(Fig.  2)  with  a  bluff  bocy  similar  to  that  described  above  centered  in 
the  strai^t  part  of  the  tube.  It  again  has  a  diameter  of  2  cm, 
althou^  a  larger  4.5  cm  bluff  body  is  also  used  for  seme  of  the 
calculations  to  determine  the  effect  of  the  bluff  boefy  size  on  the 
plasma  behavior.  For  this  configuration,  the  micrewaves  are  not  set 
in  a  standing  wave  pattern  but,  rather,  are  allcwed  to  propagate 
throu^  the  plasma  region.  Hie.  absorption  efficiencies  observed  for 
this  configuration  are  semevhat  Icwer  than  for  the  two  standing  wave 
cases.  Two  different  dewnstream  and  cenditions  are,  however,  used  for 
the  microwaves  in  this  geemetry  to  again  identic  the  important 
effects.  Hie  first  end  conditiesp  is  a  fully  transmitting  condition, 
vAiile  the  second  is  a  fully  reflecting  condition..  Hie  presence  of  a 
reflecting  end  cxmditicsi  gives  rise  to  a  weak  standing  wave  pattern 
that  in  many  ways  behaves  like  a  resonant  cavity.  Once  again,  flew  in 
this  configuration  is  from  top  to  bottem. 

2.2  Fluid  Dynamic  Model 

Hie  equations  of  moticxi  that  govern  the  flew  of  gaseous  helium 
are  the  cxnpressible  Navier-Stokes  equations.  Because  of  the  lew  Mach 
numbers  enexuntered  in  the  flew,  we  use  the  lew  Mach  number 
formulation  described  by  Choi  and  Merkle.  Hie  lew  Mach  number  version 
of  the  equations  possess^  well-ocxiditioned  eigenvalues  at  arJiitrarily 
low  Mach  numbers  and  thus  ensures  robust  cxaivergence  of  time¬ 
marching  algorithm.  For  purposes  of  the  numerical  algorithm,  the 
equations  are  written  in  the  ootple^  vficter  form  including  the 
unstea<^  terms; 
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Here,  p,  u,  v,  p  and  T  re^jresent  the  densitr^:,  x  and  y  velocity 
oonpoi>ents,  the  pressure  and  tenperature,  re^)ectiyej.y.  e  is  tte 
total  energy  per  unit  voltane  and  is  given  ky  e  =  pcyT  +  l/2p  (u*  +  v^) . 
Ihe  source  terms  are  writtai  as  tSiO-vectors^  Hi  and  H2.  Ihe  vector  Hi 
contains  the  standard  fluid  ^namic  source  t^ms  that  arise  h^xaijse  of 
the  axisynmetric  Oixirdinates.  Ihe  vector  H2  oontoins  the  source  tenos 
that  arise  fron  intoractic3ns  with  the  electronagnetic  radiation. 

These  inbcracticn  terms  are  ocnposed  of  the  IMS  values  of  the  electric 
(E)  and  magnetic  (H)  fields  and  the  ^ectxic  current  density  (J) ,  The 
energy  egua^on  contcdns  the  heat  source  frcm  electronagrietic 
radiatiOT  given  by  J  ♦  E.  In  the  mgaentum  equa^tigns,  the  source  terms 
arising  frcm  buoyancy  and  the  lorento  forces  (J  x  B)  are  included.  At 
the  present  time,  r^ation  losses  frcm  the  plasma  are  not  accounted 
for. 

The  numerical  solution  of  Ec^.  (1)  is  obtained  by  using  Euler 
ioplicit  discretization  in  time  along  with  cented.  differences  in 
space  for  both  the  inviscid  and  viscous  fliixes.  Douglas-<;unn 
cpproximate  factorization  is  us4ad  for  the  efficient  solution  of  the 
resulting  matrix  operator.  The  soluticar  procedure  involves  a  block- 
tridiagcaial  inversicp  for  each  coordinate  direction  at  eadi  time  (or 
iteration)  step.  Wtc  time-mardiing  algbritto  may  be  written  in  the 
follcwing  'delta'  form; 
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vAiere  R  is  the  steacfy  state  versicai  of  Ecji.  (^ .  Bquatioi  (2)  ensures 
that  the  steacfy  equation  is  satisfied  as  tQ  tends  te  ^ero.  The  other 
terns  in  Ec^.  (2)  are  defined  as  follows: 
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Here,  and  are  the  x  and  y  derrmtives  of  the  viscous  vector  Qy- 
^Aiich  is  define  as: 


Qv  =  [0,  u,  V,  T]t. 

In  general,  a  small  amount  of  fourth-order  artificial  ^i5fiipation  is 
added  to  prevent  odd-even  q?litting  of  the  solution. 

2.3  Maxwell  Equation  Solxitions 

The  Maxvell  equaticxis  for  TM  waves  ray  be  written  in  vector  form 
an^ogous  to  tee  fluid  (fynamic  equations:: 
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In  Eqn.  (3)  ,  the  first  three  rows  represent  the  standard  Maxv«ll's 
curl  equations  relating  the  magnetic  field  (Hig)  to  the  electric 
field  (Ejc  and  Ey) .  Ihe  last  two  equations  are  equations  for  the 
electric  curreric  density  derived  fron  momentum  consideraticxis  for  the 
conducting  electrrais  in  the  plasma.  Ihe  electrical  properties  of  the 
medium,  namely  the  electrical  oonduc±ivity  6^,  the  collision 
frequency  v^/  the  electrical  permittivitY  e  and  the  magnetic 
permeability  ,  are  calculated. 

Ihe  Maxwell  equaticais  are  solved  in  the  time  domain  ty  adopting 
an  ejq)licit  time-marching  procedure  that  is  seccsTd-rorder  accurate  in 
time  and  ^ace.  The  algorithm,  is  similar  to  that  ocnceived  by  Yee  [4] 
and  may  be  written  as  the  following  two-st^  procedure: 
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The  subscript  e  refers  to  electric  field  ccnponents  and  the  current 
densities  and  h  refers  to  magnetic  field  ccnponents.  The  first  st^ 
(Ecji.  4a)  advances  the  magnetic  field  solution  in  time  vhile  the 
second  step  vpdates  the  electric  field  solution. 

The  above  soluticxi  procedure  yields  an  unstea(fy  electrcroagnetic 
field.  Hcwever,  in  tarns  of  fluid  time  scales,  the  field  is  steady 
and  is  r^resented  by  the  BMS  values  of  the  time  varying  ccnponents. 

In  the  ccxpled  solution  procedure,  for  each  NS  (Navier-Stokes)  time 
step,  we  calculate,  the  EM  field  over  s^eral  cycles  (between  1  ai>d  10 
cycles)  and  then  tiae-average  it  to  cbtein  the  EWS  field.  Typically, 
more  cycles  were  required  per  NS  iteratioi  step  toward  the  start  of 
the  calculation.  After  about  100  iterations,  it  was  gener^ly 
sufficiart  to  ocnpute  a  single  EM  cycle  per  NS  step.  Because  of  the 
explicit  nature  of  the  Maxwell  algorithm,  the  procedure  was  found  to 
be  fairly  efficient.  Details  of  the  convergence  of  the  algorithm  and 
ccnputer  times  were  r^rted  earlier. 

2.4  Bcundarv  Conditions 

Inflow  and  outflow  boundary  conditions  for  both  the  fluid  dynamic 
and  the  Ma.i(well  equations  are  enforced  by  a  vectorized  inplementatioh 
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of  the  Method  of  Characteristics.  An  overview  of  the  procedures  is 
given  belcw. 

2.4.1  Boundary  cnTY^itinns  for  the  Navier-Stokes  Equations 

Bie  flxiid  dynamic  quantities  an  the  boundaries  are  ccoputed  by 
canbining  Method  of  Characteristic  procedures  with  physically 
realistic  boundary  oonditic*is  [5] .  At  the  inflow  boundary,  the  total 
pressure,  the  stagnaticai  tenperature  and  the  flew  angle  are  ^)ecif ied 
alcaig  with  an  tpstream  running  characteristic.  At  the  cutflcw 
boundary,  the  pressure  is  specified  along  with  three  cutrunning 
characteristics.  Synmetry  and  anti-symmetry  oonditiais  are  enforced 
ai  the  axis  of  symmetry,  vhile  no  slip  conditicxis,  an  ^prepriate 
tenperature/heat  flux  boundary  condition  (as  specified  in  the 
discussion)  and  the  normal  mcmentum  equation  are  enforced  on  the  walls. 
These  boundary  procedure  proved  to  be  very  robust  and  provided 
physically  realistic,  well-bdiaved  solutions  at  all  boundaries. 

2.4.2  Boundary  Conditions  for  the  Maxwell  Equations 

Boundary  oonditicxis  for  the  Maxwell  equations  were  likewise 
enforced  by  means  of  the  Method  of  Characteristics.  Ifexwell's 
equations  are  a  set  of  hypeibolic  equations,  and  like  any  such  set, 
they  can  be  transformed  to  characteristic  form  ty  preraultiplying  the 
equations  throu^  fcy  the  modal  matrix  of  the  eigenvectors  of  the 
equation  set.  For  the  x-direction,  we  get: 

+  A  ^  I*  =  (5) 

\diere  is  the  modal  matrix  for  the  x-direction,  is  a  diagonal 
matrix  contedning  the  eigenvalues  of  the  system  and  ^  is  the 
characteristic  vector  in  the  x-direction  that  is  defined  as 
and  is  given  by: 


/  +  yJjEy)y 


Qr  = 


Exy 


The  characteristic  vector,  Q,  contains  the  forward  and  baeikward 
traveling  waves  or  the  Riemann  invariants  of  the  system.  The 
characteristic  form  of  the  equations  is  inportant  for  speci^ing 
boundary  conditions  particularly  vAiere  the  inlet  waveguide  feeds  tl^ 
microwave  energy  into  the  cavity.  Here,  the  inocming  characteristic 
r^resents  the  pewer  that  is  input  into  the  cavity  (or  the  incident 
power)  and  has  to  be  specified  as  a  boundary  condition.  The  cutgoing 
characteristic  represents  the  pewer  that  is  reflected  back  by  the  load 
to  the  source.  The  valxie  of  this  characteristic  depends  on  the 
soluticxi  of  the  field  in  the  cavity  and  the  appropriate  characteristic 
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equaticai  is  xised  instead  of  a  boundary  condition.  If  ciiaracteristic 
variables  are  not  \ised  at  boundaries,  it  becomes  necessary  to  invoke 
an  a^dcward  iterative  procedure  to  control  the  amount  of  power  absorbed 
in  the  plasma,  and  it  is  inpossible  to  distinguish  incident  power  from 
reflected  power  so  efficiency  cannot  be  calculated. 

The  boundary  ccxiditiaTs  for  the  waveguide  soluticais  deserve 
^jecial  menticxi.  Previous  research  in  the  area  of  radio  frequency 
plasma  ^)ecified  the  tangential  ocoponent  of  the  electric  field  at  the 
boundaries.  While  this  is  the  physically  correct  boundary  condition 
for  a  reflecting  wall  (;Aiere  the  tangential  electric  field  is  zero) , 
it  is  not  proper  for  the  power  inlet  and  outlet  secticsis  of  the 
waveguide.  Ihis  is  because  the  value  of  the  electric  field  at  these 
boundaries  d^)ends  on  the  amount  of  pcwer  dissipated  in  the  plasma  and 
cannot  be  known  a  priori.  In  cin  earlier  paper,  we  shewed  that  the 
correct  quantities  to  keqp  fixed  at  these  boundaries  are  the  Riemann 
variable  (or  characteristic) .  At  the  inlet  section,  the  forward 
rtmning  characteristic  r^resents  the  incident  power  and  is  specified 
as  a  boundary  boundary  condition  vhile  the  backward  running 
characteristic  r^resents  the  reflected  pcwer  cenponent  and  is 
calculated  from  the  solution.  At  the  exit  section,  the  outgoing 
characteristic  r^resents  the  pcwer  transmitted  thixui^  the  waveguide 
\diile  the  inward  running  characteristic  represents  ary  pcwer  that  is 
incident  at  this  end.  Earlier,  we  discussed  two  kinds  of  conditions 
that  may  be  applied  at  the  exit  of  the  waveguide  —  transmitting  and 
reflecting.  For  the  waveguide  with  a  transmitting  exit  condition,  the 
inward  running  characteristic  is  zero.  For  the  waveguide  with  a 
reflecting  end  wall  at  the  exit,  the  inward  running  characteristic  is 
equal  to  the  negative  of  the  outgoing  characteristic  (i.e. ,  the 
transmitted  pcwer  is  reflected  back:  teward  the  plasma) . 

Ihe  explication  of  the  locally  one-dimensional  Method  of 
Characteristics  yields  the  follcwing  ejqjressioTS  for  the  forward  (+) 
and  backward  (-)  running  characteristics: 

9"  =  2  +  W  '  9'  =  ^  (Ey  Z„Hg)  (6) 

vAiere  Zq  is  the  characteristic  iitpedance  given  by/ 

As  a  cheeik  of  our  boundary  procedure  and  to  demonstrate  the 
differences  between  the  two  undistorted  electric  fields,  we  plot  the 
Maxwell  solutiens  for  a  lossless  waveguide  in  Fig.  3  for  both 
transmitting  and  reflecting  exit  conditions.  Figure  3a  shews  the  BMS 
value  of  the  axial  electric  field  for  the  waveguide  with  a  reflecting 
end  wall.  The  characteristic  equatiois  yield  the  obvious  ocxidition 
Ey  =  0  for  this  case.  Ihe  result  clearly  shews  a  standing  wave 
pattern  in  the  waveguide.  Since  no  power  is  dissipated  in  the 
waveguide,  ^1  the  incident  pcwer  is  transmitted  throu^  to  the  end 
wall,  vihere  it  is  reflected  back  teward  the  microwave  source.  It  is 
also  evident  that  the  exit  wall  acts  as  a  node  of  maximum  intensity. 
One-half  wavelength  ipstream  of  the  end  wall  is  a  second  node  of 
maximum  intensity.  As  mentiaied  earlier,  this  node  would  serve  as  an 
ideal  locatiai  for  the  plasma  and  good  plasma  stability  is  anticipated. 
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(We  will  show  that  this  is  indeed  the  case  in  the  Results  sectioi  of 
this  paper.)  In  addition,  the  location  of  the  waveguide  end  wall  may 
be  adjusted  in  order  to  fine  tune  the  location  of  the  plasma. 

Ihe  waveguide  with  the  transmitting  end  condition  is  shewn  in 
Fig.  3b.  Surprisingly,  this  plot  of  the  axial  electric  field  again 
shows  a  standing  wave  pattern,  edtheu^  it  is  weaker  than  the  standing 
wave  in  Fig.  3a.  Biis  means  that  not  all  of  the  incident  power  is 
transmitted  throu^  the  exit  section;  ratter,  a  percent  of  the  power 
is  reflected  at  this  boundary  (22%  in  this  case) .  Ite  reason  for  this 
is  the  use  of  Icxally  one-dimensional  Method  of  Characteristics  at 
this  (two-dimensional)  boundary.  For  the  axisymnetric  waveguide,  the 
one-dimensional  boundary  procedure  cannot  eliminate  reflection 
ocnpletely.  To  ^ecify  this  boundary  oxndition  correctly,  it  is 
necessary  to  enplcy  the  two-dimensional  Method  of  Characteristics.  In 
fact,  this  is  a  familiar  problem  in  cxmputaticml  fluid  dynamics  and  a 
formidable  one.  Fortunately,  however,  for  microwaves  traveling  in  a 
lossless  waveguide,  there  is  an  exact  solution.  Using  this  exact 
solution,  it  becomes  evident  that  the  cxsrrect  characteristic 
quantities  are  given  by: 

9'  =  i(Ey'*-W'  9'  =  2<Sy-W  <’> 

^Aiere  Zjm  is  the  ciaraoteristic  impedance  for  TM  waves  in  a  waveguide 
and  is  given  by/  n'/E  -  kVw^  Here,  u  is  the  frequency  of  the 
microwaves  and  k^  is  the  wavenunber  in  the  radial  direction. 

Figure  3c  shows  the  result  for  the  transmitting  waveguide  using 
Ekji.  (7) .  Now  the  field  pattern  is  carpletely  changed.  Ihe  standing 
wave  solution  is  no  longer  evioient;  instead,  the  field  lines  run 
parallel  to  the  waveguioie  axis  with  the  maximum  oxxurring  at  the 
centerline.  Ciis  is  the  oxjrrect  solution  for  traveling  waves  in  the 
waveguide.  The  power  fluxes  of  the  calculation  also  shew  that  all  the 
incident  power  is  transmitted  throu^  the  exit  and  none  is  reflected. 
It  should  also  be  noted  that  for  this  traveling  wave  solution  there  is 
no  preferred  axial  loxaticn  for  the  plasma;  henoe,  the  need  for  a 
stabilization  oJevice  such  as  a  bluff  boofy  is  evident. 
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3.0  RESUmS 


A  sunnary  of  parametric  results  cxi  the  three  microwave  thruster 
configurations  is  presented  in  this  section.  We  begin  with  the 
results  on  the  bli;iff-bo(^  stabilized  resonant  plasma,  thai  move  to  the 
&ee-floating  resonant  cavity  plasmas  and  finally  end  with  the 
waveguide  plasmas. 

3.1  Bluff-Body  Stabilized  Resonant  cavity  Pia.c3was 

Ihe  geonetry  for  the  configuration  corresponding  to  the 
ea^grimental  results  of  I%a±  1  is  givai  in  Fig.  lb.  A  second 
geometrical  configuration  with  a  larger  bluff  boc^  shape  and  a 
scme^tot  different  flow  geometry  is  given  in  Fig.  4.  We  begin  with 
solxitions  for  the  larger  bluff  body  (Fig.  4) .  First,  we  show  a  set  of 
results  for  cold  flow  conditions  to  fix  the  details  of  the  flowfield. 
Figure  5  shows  the  ocnvergenoe  rates  for  a  sequax^e  of  cold  flow 
calculations  in  the  geometry  of  Fig.  5  for  a  series  of  Reynolds  nunober 
conditions  that  more  than  cover  the  ^lectrum  of  the  e}^)erimental 
measurements.  Flowfield  conditions  for  Reynolds  nunbers  from  5  to 
1000  have  been  ocnputed.  As  the  convergaxoe  rates  show,  convergence 
is  good  over  all  Reynolds  nunber  levels,  with  convergence  of  three  to 
four  orders  of  magnitude  being  obtained  in  1000  iterations. 
Oorresponding  results  for  several  of  these  cases  are  given  in  Fig.  6. 
At  low  ^leeds  (low  Reynolds  nmbers)  the  cold  flow  closes  over  the 
bluff  bo^  with  no  recirculation  region,  but  as  the  ^)eed  (Reynolds 
nunber)  is  increased,  a  recirculation  region  develops  that  increases 
in  lei^.  Ihis  flowfield  is  representative  of  the  undisturbed  flow 
into  vbich  the  microwave  plasma  is  introduced.  Of  course,  the  plasma 
heating  alters  the  flowfi^d  dramatically. 

Figure  7  shows  the  convergaice  of  the  conbined  ttaxwell/N^er- 
Stoikes  equations  for  a  representative  case.  As  can  be  seal, 
convergaice  renains  good  yihen  the  plasma  is  presait.  We  do,  however, 
note  that  the  process  of  finding  an  aco^itcble  initial  condition  can 
be  challaiging,  and  the  first  solution  in  a  sequenoe  is  generally 
found  by  a  cut-and-try  method.  If  a  "hot"  zone  is  not  included  in  the 
initial  condition,  there  will  be  no  absorption.  Just  as  in  the 
experimaital  case,  the  initial  condition  mist  start  with  seme  "break¬ 
down"  region  \iMch  mist  then  be  sustained  to  obtadn  power  absorption. 


The  flowfield  solution  for  a  representative  case  is  ^xwn  in  Fig. 
8.  As  can  be  seen,  the  plasma  sits  in  the  wake  of  the  bluff  body,  but 
it  causes  a  dramatic  wealoening  in  the  size  and  strength  of  the 
recirculation  zone,  eventually  leading  to  the  ocoplete  removal  of 
recirculation  as  is  shown  later.  These  calculations  are  for  the  THOU 
mode,  and  the  plasma  forms  at  the  lower  node.  The  bluff  boc^  has  been 
situated  so  that  the  wake  region  corresponds  with  this  node. 

The  temperature  contours  for  a  series  of  differoit  microwave 
powers  are  shown  on  Fig.  9.  As  the  power  is  increased  the  plaana 
increases  in  length,  but  the  radial  growth  is  restricted  by  the  size 
of  the  bluff  boc^.  The  experimartal  measurenmits  have  shown  similar 
changes  in  the  plasma  shape  with  increases  in  the  power.  We  also  note 
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that  increased  heat  addition  causes  a  decrease  in  the  recirculation 
region,  eventually  almost  eliminating  it  at  the  higher  (4  kW)  pcwer 
level. 

Similar  tenperature  contour  plots  showing  the  effect  of  pressure 
level  are  given  on  Fig.  10.  Increases  in  pressure  caxise  diffusion  to 
becone  smaller  so  that  the  recirculation  region  beccces  stronger. 
Increased  pressures  also  lead  to  stronger  effects  of  buoyancy  vhich 
act  in  a  direction  to  augment  this.  Ihe  plasma  also  becomes  notably 
more  compact  at  the  hi^ier  pressures,  again  -in  agreement  with 
qualitative  cbservatiois  obtained  from  the  e}^}eriments. 

Representative  calculations  for  the  ejqjerimental  geometry  of  Fig. 
lb  are  presented  in  Fig.  11.  For  this  case,  the  flow  is  accelerated 
throu^  an  orifice  in  the  bottom  of  the  ohantoer  duplicating  the 
ejperimental  setup.  R^resentative  tanperature  contours  and 
streamlines  are  given  cai  Fig.  11. .  Note  that  in  this  case,  there  is  no 
recirculation  in  the  base  region,  althcu^  there  is  a  weak 
recirculation  cell  in  the  base  of  the  bulb  region  vhere  the-  narrow 
orifice  creates  a  "dead-water"  region.  The  peak  temperatures  for  this 
case  are  about  11,500  K,  a  valxoe  that  is  in  good  agreement  with 
experiment.  Again,  the  plasma  size  is  restricted  by  -the  size  of  the 
bluff  body,  so  -the  present  plasma  is  quite  narrow  in  size.  larger 
bluff  bocfy  sizes  (and  the  increased  plasm  diameters)  enable  more 
microwave  energy  to  be  added  to  the  flow  so  -that  the  average 
■tenperature  of  -the  heated  gas  increases,  "thus  providing  hi^ier 
propulsive  performance.  The  plasma  she^  noted  here  is  si^ar  to 
that  observed  in  the  experiments.  Quantitative  conparisons  of  the 
shapes  is  planned,  but  the  data  from  the  euperiments  has  not  been 
reduced  as  yet. 

A  ccnparison  between  the  predictions  and  the  esperiments  for  the 
bluff-4xxfy  stabilized  plasma  is  given  in  Figs.  12  and  13.  Figure  12 
shews  the  variation  of  the  cotpling  efficiency  with  changes  in  the 
L3SS  flow  rate  at  two  power  levels,  1  and  2  M?.  Cotpling  efficiency 
is  here  defined  as  -the  ratio  of  -the  power  absorbed  -to  -the  incident 
power.  Losses  include  only  the  reflected  pcxner,  a  quantity  that  is 
small  because  of  the  uSe  of  a  cavity  standing  wave.  Ihe  predictions 
are  in  exscellent  agreement  with  -the  measurements  for  both  pewer  levels. 
Similar  ccnpariscjns  of  the  peak  tenperatures  are  given  in  Fig.  13. 

The  ccnparison  is  again  very  good  and,  further,  we  note  that  the 
smaller  bluff  boc^  provides  nearly  a  2000  K  -tenperature  rise  as 
extpared  to  -the  larger  one.  Ihe  subject  ejperimental  measurements 
were  considered  questionable  because  -they  were  so  hi^,  especially  as 
cetpared  -to  the  free-floating  plasma,  tut  the  present  calculations 
show  that  a  hi^ier  tenperature  is  justified  because  of  -the  size  of  -the 
bluff  body.  Ihe  larger  bluff  boefy  provides  peak  tenperatures  that  are 
more  in  line  with  previously  observed  microwave  plasmas,  but  it  also 
shows  that  it  is  possible  -to  obtain  control  over  -the  pekk  tenperature. 

3.2  Results  for  the  Free-Floating  Remnant  Gavitv  Plasma 

Ihe  geometrical  configuration  for  the  free-floating  resonant 
cavity  plasma  is  given  in  Fig.  la.  Here  we  show  the  cavity,  the 
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s^jhfire-cyliixler  flow  tube,  and  the  undistorted  micrcwave  field  in  the 
1M012  mode.  Ihe  initial  results  for  the  free-floating  plasma  (Fig. 

14)  Shaw  the  placement  of  the  plasma  occurs  in  the  center  of  the 
sphere  at  the  middle  node  of  the  landisturbed  E-M  field  in  agreeanent 
with  ei^jerimenta''  ibservaticns.  Ihe  predictiois  indicate  that  the 
conditions  choGa.  the  experiments  lie  near  a  boundary  vhich 
differentiates  between  the  existence  of  and  the  absence  of  a 
recirculating  region  in  the  ^here,  but  that,  in  general,  the 
ei^jeriments  are  run  just  above  this  boundary  so  that  recirculation  is 
not  observed.  Ihe  field  lines  for  the  absorbing  electric  field  are 
given  on  Fig.  15.  Note  by  ocnpariscn  with  Fig.  la,  the  substantial 
amount  of  distorticxi  that  is  observed.  Clearly,  the  prediction  of  the 
plasna  requires  a  ooi?)led  Maxwell/Naviei>-Stoikes  sol\ition.  Ihe 
distorted  electric  field  also  gives  clear  evidence  of  the 
"skin  effect"  that  concentrates  the  field  at  the  surface  of  the  plasma 
and  leaves  the  interior  of  the  plasma  nearly  constant. 

The  effect  of  microwave  pcwer  on  the  sh^  and  location  of  a 
free-floating  plasma  in  a  strai^t  cylindrical  duct  is  given  on  Fig. 
16.  For  a  series  of  input  powers  starting  at  lew  powers,  the  plasma 
forms  at  the  bottom  of  the  three  nodes  in  the  'IM012  mode,  and 
ccxitinues  to  grow  in  size  as  the  power  is  increased  from  1.5  to  2  kW. 
At  the  3  kW,  the  plasma  begins  to  e^qpand  above  the  lowest  node  and  to 
absorb  a  small  fraction  of  the  radiation  in  the  middle  node,  and  at  4 
kW  it  juirps  to  the  top  and  middle  nodes.  Using  this  4  kW  solution  as 
the  initial  cxndition  for  first  a  3  ,kW  and  then  a  2  kW  calculation, 
the  power  in  the  plasma  is  absorbed  from  the  middle  node  of  the 
standing  wave.  These  results  show  that  the  location  of  the  plasma  in 
a  strai^t  tube  d^iends  i^xon  the  initial  condition,  i.e.,  it  d^)ends 
upon  the  manner  in  vhich  the  plasma  was  initiated.  No  ejqperimental 
results  are  available  for  a  strai^t  cylindrical  tube,  so  no 
verification  of  this  multiplicity  of  plasma  locations  has  been  shown, 
but  the  expansion  of  the  plasma  at  hi^ier  powers  to  cover  two  nodes 
has  been  observed  ej^jerimentcilly.  Similar  calculations  for  the 
qphere-cylinder  tube  in  Fig.  17  show  the  plasma  forms  at  the  center 
node  and  expands  in  size  with  increasing  iiput  power  in  agreement  with 
ejqjeriment.  These  results  vhich  are  chosen  to  match  the  e)5)erimental 
results  described  in  Part  I,  show  that  the  flowfield  in  the  spherical 
region  contains  a  large  recirculation  region  at  the  lower  power 
levels,  but  that  this  recirculation  is  weakened  as  the  power  is 
increased,  and  eventually  is  eliminated  at  the  4  kW  case.  These 
sphere-cylinder  results  also  show  the  plasma  size  expands  to  form  a 
double  plasma  at  the  hi^ier  powers,  in  agreement  with  the  experiments. 

The  effects  of  increasing  the  inlet  gas  velcwity  from  0.1  to  1 
Tty's  are  shown  in  Figs.  18  and  19  for  the  sphere-cylinder  geometry  and 
the  strcii^t  cylindrical  geemetry,  respectively.  In  both  cases,  the 
hi^rer  inlet  velocities  displace  the  plasma  downward  away  from  the 
incoming  flow.  The  hi^ier  inlet  velocities  reduce  the  ipward 
diffusion  of  the  heat  and  indicate  the  location  of  the  energy 
absorption  regions  more  directly.  For  both  geometries,  the  effect  of 
bucyaixy  is  more  pronounced  at  the  lower  velocities  vhere 
recirculation  is  observed,  while  at  the  hi^ier  velocities,  the 
recirculation  region  is  no  longer  present.  Corresponding  effects  of 
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changing  the  pressure  on  both  the  j^here^lihder  and  the  strai^t 
cylinder  are  shown  in  Figs.  20  and  21.  As  for  the  bluff-bocfy 
stabilized  plasm,  the  larger  plasms  reduce  the  size  of  the  plasm 
making  it  more  ccnpact.  Ei^erimental  measurements  have  also  noted 
this  size  change  with  increases  in  the  pressure  level.  Ihe  increased 
pressure  cilso  increases  the  size  of  the  recirculation  region  because 
the  increased  pressure  decreases  the  rate  of  diffusicsi  caiasing 
diffusive  effects  to  beocrae  more  significant  relative  to  buoyancy. 

Seme  direct  cenparisons  with  experiment  are  given  in  Figs.  22  and 
23.  Figure  22  shows  the  predicted  and  measured  peak  tenperatures  as  a 
function  of  the  pressure.  Again,  the  peak  tenperature  is  nearly 
independent  of  mass  flow  for  both  the  predictions  and  the  measurements. 
The  predictions  are  some  1000  K  Icwer  than  the  measurements,  but  the 
difference  between  the  present  free-floating  plasm  and  the  previc^ly 
discussed  bluff-bocfy  stabilized  plasm  is  well  predicted  by  the 
analysis.  The  predicted  cotpling  efficiencies  for  the  free-floating 
plasm  are  given  in  Fig.  23,  and  shew  that  the  free-floating  plasma  is 
capable  of  cocpling  almost  all  of  the  incident  energy  into  the  plasm. 
There  is  seme  fall-off  in  efficiency  predicted  at  the  lowest  and 
hi^er  pressures.  The  reason  for  this  is  that  the  present  cavity 
incorporates  a  sliding  short  that  can  be  used  to  re-tune  the  cavily 
length  in  the  presence  of  losses  (absorption) .  The  epperimentcil 
measurements  of  efficiency  are  considerably  lower  than  the  predictions 
because  th^  include  heat  losses  to  the  wall  and  electrcxnagnetic 
losses  through  the  wall.  The  experiments  were,  hewever,  conducted  in 
such  a  way  that  the  sliding  short  was  adjusted  in  the  presence  of  the 
plasm  until  the  reflected  power  was  zero,  or  nearly  so,  indicating 
the  true  coupling  efficiency  was  nearly  1000%,.  in  agreement  with  the 
predictions. 

3.3  Results  for  the  Bluff-Body  Rtab-i  1  jzed  Plasm  in  a  Wavecaiide 

3.3.1  Representative  Wavecaiide  Soliition 

To  provide  a  physical  picture  of  waveguide  heated  plasmas,  a 
typical  solution  for  the  waveguide  with  a  transmitting  boundary  is 
shown  in  Fig.  24.  The  flowfield  solution  is  given  in  Fig.  24a  and  the 
electric  field  solxition  is  in  Fig.  24b.  For  this  case,  the  inlet  gas 
velcpcity  is  3  Vs  and  the  gas  pressure  is  1  atm;  the  corresponding 
mass  flow  rate  is  1.1  g/s.  The  incident  micrewave  pewer  is  3.5  kW, 
vhile  the  diameter  of  the  bluff  boefy  used  to  hold  the  plasm  is  4.5 
cm. 


The  temperature  contours  in  Fig.  24a  show  that  the  plasm 
discharge  is  maintained  in  the  wake  of  the  bluff  boc^.  As  shewn  in 
Fig.  3,  the  electric  field  in  the  transmitting  waveguide  does  not 
offer  any  preferred  locaticxi  for  the  plasm  so  the  Icxatiai  of  the 
plasm  is  determined  purely  by  the  (fynamics  of  the  flew.  The  fluid 
(fynamics  tends  to  push  the  plasm  dewnstream  \Mle  the  micrewave  power 
pulls  it  vpstream.  The  plasm  is  stable  in  a  region  vhere  these  two 
forces  balance  out.  The  plasm  is  long  and  slender  with  its  outer 
diameter  being  defined  by  the  edge  of  the  bluff  body.  This  is  because 
of  hi^  local  gas  velcxsities  at  the  outer  edge  of  the  bluff  boeV*  Bie 
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peak  plasma  temperature  axurs  on  the  oentxal  axis  and  is  about  9150  K 
vMch  is  typical  of  helium  plasmas  of  this  size.  Ihe  streamline 
omitours  show  evidence  of  a  recirculating  region  in  the  wake  of  the 
bluff  bodn,  however,  the  recirculating  region  is  quite  weak  relative 
to  its  cold  flow  counterpart  under  similar  conditions.  Ihis  is 
because  the  strong  heat  addition  in  the  plasma  tends  to  suppress  the 
recirculation  zoie.  In  fact,  in  sane  cases,  the  recirculation  is 
totally  eliminated. 

Ihe  axial  and  radial  electric  field  contours  are  shewn  in  Fig. 

24b.  Ihe  electric  field  lines  are  seen  to  be  strongly  distorted  in 
the  presence  of  the  plasma  (oenpare  with  Fig.  3) .  Because  of  the  skin 
effect,  the  absorption  of  the  energy  occurs  at  the  surface  (or  skin) 
of  the  plasma.  Consequently,  strong  ^ectric  field  gradients  exist  on 
the  surface  vhile  the  electric  field  is  zero  inside  the  plasma.  Of 
the  3.5  kW  that  is  incident  on  the  plasma,  2.4  is  absorbed  by  the 
plasma,  1  bW  is  reflected  back  toward  the  inlet  and  0.1  ]<W  is 
transmitted  through  the  waveguide  exit.  Ihe  ootpling  efficiency, 
vAiich  is  th2  ratio  of  absorbed  to  incident  power,  is  70%  for  this  case. 
In  Fig.  24b,  the  power  reflected  by  the  plasna  creates  a  standing  wave 
pattern  ipstream  of  the  plasma. .  Dewnstream  of  the  plasma,  the 
electric  field  lines  run  parallel  to  the  waveguide  axis,  r^resenting 
the  transmitted  component.  For  the  transmitting  waveguide 
oonfiguraticai,  both  the  reflected  and  transmitted  powers  are  lost  from 
the  system. 


3.3.2  Oemoarison  of  Transmitting  and  Reflecting  Waveguides 

In  a  reflecting  waveguide,  the  power  transmitted  to  the  exit 
section  is  reflected  back  toward  the  plasma  and  thus  gets  a  second 
chance  at  being  absorbed.  It  is  interesting  to  check  if  this 
alternate  configuration  gives  better  energy  coupling  than  the 
transmitting  waveguide.  Figure  25  shows  the  results  for  a  reflecting 
waveguide  for  the  same  oonditicxis  used  in  Fig.  24. 

Ihe  reflecting  waveguide  solution  in  Fig.  25  is  remarkably 
similar  to  the  transmitting  waveguide  soluticn.  Ihe  peak  temperature 
is  now  9310  K  vhile  the  coupling  efficiency  is  68%.  Ihe  streamline, 
not  shown,  are  also  qualitatively  similar.  Ihe  electric  field 
soluticai  now  shows  evidence  of  a  standing  wave  pattern  dewnstream  of 
the  plasma  as  well  (indicating  power  refleotion  at  the  exit) . 

However,  the  reflecting  ccaiciLtion  has  made  very  little  impact  on  the 
nature  of  the  plasma.  One  reason  for  this  may  be  that  only  a  small 
amount  of  power  is  transmitted  throu^  the  plasma  in  this  case. 

Figure  26  compares  the  performanoe  of  the  reflecting  and 
transmitting  waveguides  over  a  wide  range  of  powers;  the  reflecting 
waveguide  results  are  shown  cn  the  top  (Figs.  26a  to  26d)  vMle  the 
transmitting  waveguide  results  are  shown  on  the  bottom  (Figs.  26e  and 
26f) .  Ihe  flow  conditiens  are  once  again  the  same  as  before  (1  atm 
and  3  Bv/s  inlet  velocity) .  For  the  reflecting  waveguide  (Figs.  26a  to 
26d) ,  the  plasma  is  maintciined  stably  in  the  wake  of  the  bluff  body 
over  the  range  of  powers.  As  the  power  is.  reduced  from  5  kW  to  1  kW, 
the  plasma  ^ihks  in  size.  Ihis  effect  is  more  noticeable  in  the 
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pla^  length  than  in  its  radius  because  the  bliiff  body  restricts  the 
reuiial  size  6f  the  plasma.  Hie  peak  plasma  tenperature,  hcwev^, 
stays  constant  at  about  9300  K.  Of  particular  interest  is  the  double 
plasma  -Uiat  is  starting  to  form  at  the  5  kW  condition.  Such  double 
plasmas  have  also  been  observed  experimentally.  When  the  power  is 
reduced  below  1  W/i,  the  plasma  is  eventually  extinguished.  Such  a 
pc^^  threshold  exists  because,  at  low  powers,  the  electrcais  start 
reccmbining  more  rspidly  than  they  are  created  by  ionizaticxi.  In  this 
case,  this  minimum  pcwer  was  about  800  W,  \»hich  is  similar  to 
eiqieriiiiental  values. 

An  interesting  observation  for  the  reflecting  waveguide  is  that 
the  cotpling  efficiency  of  the  plasma  increases  as  the  power  is 
reduced,  starting  fran  58%  (at  5  kW)  to  95%  (at  1  kW) .  This  is 
probably  because  at  hii^  powers  the  plasma  elongates  in  size  but  is 
restrained  frcm  growing  radially  by  the  bluff  bocfy.  As  a  result,  the 
plasma  reflects  more  at  hi^ier  powers  and  absorbs  less.  Hiis  result 
is  quite  different  from  the  results  for  free-floating  resonant  cavity 
plasmas  vhere  coipling  efficiencies  vp  to  100%  are  obtainable  at  all 
pcwer  levels. 

Hie  transmitting  waveguide  results  for  similar  conditions  are 
shewn  in  Figs.  26e  and  26f.  At  the  hi^er  powers  (5  kw  and  3.5  kW) , 
the  plasma  discharge  is  qualitatively  similar  to  the  reflecting 
waveguide  plasma.  But  at  powers  lower  than  about  3  kW,  it  was 
iirpossible  to  maintain  a  stable  plasma.  Hiis.is  because  at  these 
lower  powers,  the  convection  d/namics  of  the  flow  deminates  over  the 
microwave  pewer  and  blows  the  plasma  away  frcm  the  wake  of  the  bluff 
bocty.  In  other  words,  the  mass  flow  rate  associated  with  the  plasma 
is  higher  than  the  plasma  propagation  ^)eed  at  this  ccxidition.  Hius, 
the  transmitting  waveguide  can  maintain  stable  plasmas  only  above  this 
power  threshold  of  3  kW.  Of  course,  for  lower  gas  flow  rates,  this 
threshold  would  be  lowered;  but,  for  the  flow  rates  of  interest  here, 
no  significant  change  in  the  threshold  power  was  observed.  Haerefore, 
the  reflecting  waveguide  configuration  can  maintain  stable  plasmas 
over  a  wider  range  of  powers  than  the  transmitting  waveguide. 

Hie  inproved  performance  of  the  reflecting  waveguide  arises 
becaose  it  absorbs  energy  in  the  standing  wave  mode  \Aule  the 
transmitting  waveguide  absorbs  in  the  propagating  wave  mode.  For  the 
propagating  jaode,  the  electric  field  contours  do  not  provide  a 
preferred  '  .  aticai  for  the  plasma,  vtile,  for  the  standing  wave  mode, 
the  nodes  oi  iuaximum  field  intensity  do  provide  such  a  location.  For 
the  reflecting  end  wall,  one  such  noefe  occurs  one-half  wavelength 
vpstream  of  the  end  wall.  Hie  octtoined  presence  of  this  node  and  the 
wake  of  the  bluff  body  ccxTtritwtes  to  increased  stability  of  the 
plasma  over  a  wide  range  of  irput  powers.  In  fact,  experiments  have 
indicated  that  it  is  senetimes  possible  to  maintain  the  plasma  even  in 
the  absence  of  the  bluff  body.  For  the  parametric  studies  in  the  rest 
of  this  r^rt,  we  ooncxiiitrate  on  the  reflecting  configuration. 
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3.3.3  Locaticai  of  thg  Pla«gna 

Ihe  locatioi  of  the  plasma  in  the  reflec±ing  waveguide  is 
influenced  ty  the  locatioi  of  the  maxinun  node.  Ihus,  ty  adjusting 
the  locaticn  of  the  reflecting  end  wall,  it  should  be  possible  to 
control  the  precise  location  of  the  plasma.  Figure  27  shews  the 
results  when  the  reflecting  end  wall  is  positdwied  at  two  different 
dcwnstu:eam  locations  —  in  Fig.  27a,  the  distance  between  the  base  of 
the  bluff  bo(fy  and  the  end  wall  is  3  (vAiere  is  the  bluff  boc^ 
diameter) ;  \iAiile,  in  Fig.  27b,  this  separation  distance  is  4  E^.  In 
Fig.  27a,  the  plasma  is  pushed  vp  against  the  bluff  bo<fy,  \Aiile  in 
Fig.  27b,  the  plasma  is  detached  freo  the  bluff  bo(fy.  Clearly,  this 
movement  of  the  plasma  is  caused  by  the  migration  of  the  maxiimm  node 
as  the  end  wall  is  moved  further  dewnstream. 

Ihe  streamline  contours  for  the  two  conditions  are  also  shown  in 
Fig.  27.  The  recirculation  region  is  seen  to  grow  as  the  plasma  is 
moved  dewnstream  of  the  bluff  boc^.  Ihe  heat  transfer  to  the  bluff 
boefy  is  influenced  by  the  strength  of  this  recirculating  region.  In 
Fig.  27a,  the  heat  conducted  from  the  plasma  to  the  bluff  body  is 
quite  hi^  (0.46  oorresponding  to  about  15%  of  the  incident  power) . 
In  Fig.  27b,  the  heat  loss  to  the  blviff  bocty  is  only  0.04  kW  (or  1.3% 
of  the  incident  power) .  The  condition  in  Fig.  27a  is  clearly 
vmdesirable.  The  hi^  heat  loss  to  the  bluff  boc^  results  in  poor 
energy  conservation  and  low  overall  thermal  efficiency.  Furthermore, 
hi^  heat  load  may  also  result  in  degradation  of  the  matericil  of  the 
bluff  body.  Thus,  the  control  of  the  plasma  location  that  is  made 
possible  with  the  reflecting  waveguide  can  ensure  better  overall 
perfontance. 


3.3.4  Effect  of  Bluff  Body  Size 

We  observed  in  Fig.  24  that  the  radius  of  the  plasma  was  defined 
by  the  size  of  the  bluff  body.  It  is  interesting  then  to  see  how  the 
plasma  characteristics  change  \4ien  a  smaller  blxiff  body  is  enplcyed. 
Figure  28  shows  the  teirperature  centours  and  streamlines  for  a  plasma 
(incident  power  =  3  W?,  1  atm  pressure,  3  u/s  inlet  velocity)  that  is 
stabilized  by  a  2  cm  diameter  bluff  bo^.  Ihe  peak  plasma  tenperature 
is  now  significantly  higher  at  about  11,600  K.  Biis  hitcher 
tenperature  is  a  consequence  of  the  smaller  size  of  the  plasma.  Such 
an  effect  has  been  previously  observed  both  experimentally  and 
theoretically.  Hiis  tenperature  is  also  in  agreement  with  earlier 
measurements  of  waveguide  plasnas,  ^hich  enployed  a  bluff  boefy  of 
similar  size.  Fmdhermore,  the  tenperature  contours  show  the  plasma 
very  closely  attached  to  the  base  of  the  bluff  bocy.  The  streamline 
contours  show  no  eviderxxe  of  recirculation  in  the  wake  indicating  that 
the  high  heat  addition  rate  in  this  region  has  carpletely  suppressed 
it. 


The  coupling  efficiency  for  this  calculation  is  47%  cenpared  to 
about  68%  for  the  4.5  cm  bluff  body.  It  is  evident  that  the  size  of 
the  plasma  restricts  its  absorption  Cc^^city  and  causes  a  large 
portion  of  the  incident  power  to  be  reflected  back.  As  in  the  earlier 
case,  the  size  of  the  bluff  body  controls  the  radial  exteiit  of  the 
plasma  and,  therefore,  plays  a  crucial  role  in  determining  the 
efficierxy  of  energy  coupling.  In  addition,  there  is  also 
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(xnsidexable  heat  loss  to  the  bluff  body  because  of  the  close 
proximity  of  tl^  plasma  to  the  bluff  botfy  (6.22  kW  or  7%  of  the 
incident  power)  resulting  in  hi^  heat  loa^  to  the  bli:iff  bod^  as  v^l. 
This  heating  o^d  again  be  oontroll^  by  fixing  the  locaticn  of  the 
bluff  boc^  wit^  re^sect  to  the  standing  weive  node. 

Vfe  have  established  sane  of  the  basic  characteristics  of  plasmas 
heated  in  a  waveguide  ocnfiguraticn.  In  tiie  remainder  of  the  r^rt, 
we  carry  cut  detcdled  parametric  studies  to  establish  hew  the  plasma 
characteristics  change  under  different  c^^erating  conditions.  The 
plasma  characteristics  examined  .are  the  size  of  the  plasma  (length  and 
width) ,  the  size  of  the  recirculatiai  regicxi,  the  coupling  and  thermal 
efficiencies  of  the  plasma  and,  finally,  the  peak  and  average  gas 
tenperatures.  The  results,  sumoairized  in  Figs.  29  to  33,  are  plotted 
against  the  mass  flow  rate  of  helium.  Mass  flow  rate  may  be  varied 
two  ways  —  fcy  changing  the  inlet  gas  velocity  and  by  changing  the  gas 
pressure.  We  consider  both  these  effects  ind^jendently.  Further,  the 
results  are  presented  for  two  incident  power  levels  (2  kW  and  4  kW) . 
Finally,  in  additioi  to  the  results  for  the  4.5  an  bluff  boc^,  sane 
results  are  shewn  for  the  smaller  bluff  bex^  size  (2  cm  diameter)  as 
well. 


Similar  parametric  studies  are  curroitly  being  performed 
e}^)erimentally  with  the  epparatus  shown  in  Fig.  2.  At  the  present 
time,  very  little  consolidated  data  are  available  for  direct 
conparison.  Earlier  ejperiments  for.  a  smaller  flew  geemetry  and  a 
rectangular  waveguide  have  yielded  plasma  tenperatures  of  12,000  K  and 
coupling  efficiencies  of  60  to  80%;  these  data  are  in  qualitative 
agreement  with  the  present  results.  In  addition,  many  of  the  trends 
predicted  ccnputationally  have  been  corroborated  experimentally. 

3.3.5  Size  of  the 

The  size  of  the  plasma  is,  naturally,  of  importance.  We  have 
alreacfy  seen  that  the  size  of  the  plasma  influences  the  coupling 
efficiency  of  the  plasma.  Rirthernore,  the  heat  loads  on  the  walls 
and  cx\  the  bluff  boc^  also  d^end  on  plasma  size.  Figure  29  shews  the 
plasma  length  and  the  mean  plasma  diameter  plotted  against  mass  flew 
rate.  Both  quantities  are  normalized  with  reject  to  the  size  of  the 
bluff  boi^.  Four  sets  of  data  are  shewn.  Here,  as  well  as  in  Figs. 

30  to  33,  the  circles  r^resent  the  case  vhere  the  mass  flew  rate  is 
varied  by  changing  the  gas  velocity  (2  kW  pewer,  1  atm  pressure) ;  the 
triangles  are  for  a  hi^ier  power  of  4  kW,  1  atm  pressure  and  different 
velocities;  the  squares  are  again  for  the  Icwer  pewer  of  2  kW  with  the 
gas  velocity  maintained  at  0.75  ny's  iihile  gas  pressure  is  varied;  and 
finely,  the  stars  shew  the  calculations  for  the  2  an  bluff  boefy  for 
two  gas  velocities  (1  atm  pressure,  3  kW  pewer) . 

The  most  noticeable  effect  in  Fig.  29  is  the  increase  in  the  size 
of  the  plasma  with  power.  The  4  kW  plasmas  are  seen  to  be 
significantly  longer  (by  about  50%)  and  semevhat  wider  (by  about  20%) 
than  the  2  kW  plasmas.  As  noted  earlier,  the  restriction  in  radieil 
growth  is  due  to  the  presence  of  the  bluff  boefy.  Increasing  the  gas 
velocity  reduces  the  size  of  the  plasma  slightly  in  roost  cases  (excqpt 
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the  length  of  the  4  kW  plasna  vMch  giX3ws  sooeiitot) .  This  is  a 
Reynolds  number  effect  —  hi^ier  gas  velocities  mean  high^  Reynolds 
number  and,  hence,  less  diffusion  of  thermal  energy.  The  plasma  thus 
becxtnes  more  ocnpact  at  hi^xr  velociti.es.  Increasing  pressure  also 
increases  the  R^nolds  number  and  a  sli^t  reduction  in  plasma  size  is 
evident.  The  results  for  the  smaller  bluff  bo(fy  show  similar  plasma 
sizes  with  re^»ct  to  the  smaller  diameter.  In  terms  of  actual  size, 
these  plasmas  are  thinner  and  about  as  long  ^  the  4.5  cm  bluff 
body  plasmas. 

3.3.6  Size  of  the  ppr.iTxjjlation  Recrion 

Figure  30  shews  the  variation  in  the  length  of  the  recirculating 
zone  in  the  wake  of  the  bluff  boefy.  As  mentioned  earlier,  the 
recirculating  zone  influ^xes  the  rate  of  heat  transfer  to  the  bluff 
body.  Depending  on  the  exact  location  of  the  plasma  with  re^)ect  to 
the  bluff  body,  the  size  of  the  recirculation  zcx>e  varies  and  so  does 
the  heat  transferred  to  the  bl^f  body.  When  the  plasma  is  in  close 
proximity  to  the  base  of  the  bo(ty,  the  recirculation  regicai  is  totally 
suppressed  (as  in  Fig.  28)  and  there  is  a  significant  heat  load  cn  the 
bluff  bocty.  Figure  30  sh^  that  the  recirculaticxi  regiai  is 
suppressed  at  low  flow  velocities  and  low  pressures,  that  is,  at  low 
mass  flow  rates.  Increasing  the  mass  flow  rate  by  increasing  gas 
velocity  or  pressure  elongates  the  recirculating  region  (the  effect  of 
pressure  is  sli^tly  greater) .  This  is  aga^  a  Reynolds  number  effect. 
At  lew  Reynolds  numbers,  diffusion  of  thermal  energy  is  significant 
and  the  tenp^ture  contciurs  extend  over  a  larger  regie*!  around  the 
plasma;  the  larger  this  region,  the  more  recirculation  is  suppressed. 
Also,  Fig.  30  shows  that  at  hitler  powers  the  recirculating  regioi  is 
suppressed  over  a  bigger  range  of  mass  flows  (or  Reynolds  nunbers) . 
Clearly,  at  hi^ier  powers  the  intensity  of  heat  addition  is  great^ 
and  the  plasma  is  bigger,  both  of  vdiioh  result  in  a  stronger 
suppiessic*!  of  recirculation. 

3.3.7  Oouplincr  Efficiency 

The  ocxp>ling  efficiency  of  the  plasma  is  defined  as  Idle  percent 
of  the  incident  power  that  is  absorbed  by  the  plasma.  Figure  31  shews 
the  coupling  efficiency  plotted  against  mass  flow  rate.  There 
eppears  to  be  very  little  d^jendenoe  on  the  £Low  rate.  For  the  2  kW 
plasma,  ooupling  efficiencies  are  about  80%  \^!ile  for  the  4  kW  plaana, 
they  are  about  65  to  70%.  The  lower  coupling  efficiencies  at  high^ 
powers  was  noted  earlier  and  may  be  attributed  to  the  restriction  in 
the  radial  growth  of  the  plasma  by  the  bluff  boefy.  This  cilso  eipledns 
vhy  the  coupling  efficiencies  for  the  2  can  bluff  body  plasmas  (3  kW) 
are  less  than  50%. 

3.3.8  Overall  Thermal  Efficiency 

Once  the  plasma  absorbs  the  radiaticai,  it  loses  energy  to  the 
surrounding  by  conc^iction  and  radiation.  It  is,  therefore,  ispertant 
to  check  hew  much  of  the  energy  is  actually  retoined  by  the  gas.  'Ihe 
overall  thermal  efficiency  is  defined  as  the  total  eiiergy  in  the  gas 
at  the  exit  section  of  the  flow  tube  divided  ty  the  incident  power. 
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By  definition,  it  \dll  be  less  them  the  cxwpling  efficiency  by  exactly 
the  amount  of  energy  that  is  lost  to  the  side  wills  and  the  bluff  boc^. 
Figure  32  presents overall  thermal  efficiency  versus  mass  flow 
rate.  Increasing  ness  flow  rate  (either  by  increasing  pressure  or 
flow  velocity)  increases  Reynolds  number,  decreases  therma],  diffusion, 
reduces  plasma  size,  increases  the  size  of  the  recirculatioi  region 
and  therdy  increases  the  thermal  efficiency.  The  overall  efficiency 
is  also  seen  to  be  hi^ier  for  the  2  bW  plasma  than  for  the  4  ysf  plasma. 
This  is  clearly  because  of  the  larger  size  of  the  4  ]?W  plaana.  Figure 
32  also  shows  that,  at  hi^  encu^  flow  rates,  it  is  possible  to 
achieve  thermal  efficiencies  i?)  to  75%,  \Auch  is  nearly  equal  to  the 
ootpling  efficiency  for  those  cases.  Iftider  these  conditicais, 
practically  all  diffusion  losses  are  eliminated.  Ihere  is,  however,  a 
price  to  pay  for  this  as  will  be  evident  in  the  next  section. 

3.3.9  Gas  Temperatures 


It  is  irportarit  to  ensure  that  the  incident  radiation  is  absorbed 
efficiently  and  retained  by  the  gas.  It  is  also  iirportant  to  check 
how  nudi  the  gas  is  heated  by  the  radiatiai.  It  is  the  latter  factor 
that  determines  the  performance  of  the  thruster  device  (the  hi^ier  the 
gas  temperature,  the  hic^ier  the  thrust  and  ^)ecific  impolse  of  the 
thruster) .  Figure  33  presents  both  the  peak  gas  temperature  that 
occurs  in  the  central  core  of  the  plasma  and  the  average  gas 
temperature  at  the  exit  of  the  flow  tube.  The  peak  gas  temperature  is 
a  useful  parameter  for  characterizing  plasmas,  especially  since  this 
is  vhat  is  measured  most  frequently  in  experiments.  Ihe  average 
temperature  is,  however,  more  important  since  it  determines  the 
performance  of  the  thruster.  Figure  33  shews  that  pe^  temperatures 
are  around  9000  to  10,000  K  for  the  4.5  cm  bluff  boefy  and  about  11,500 
K  for  the  2.  cm  bluff  bo^.  The  peak  temperature  is  relatively 
insensitive  to  flew  conditicais;  but  a  sli^t  increase  is  noticed  with 
increasing  gas  flew  rate.  Ihis  latter  effect  may  be  ejplained  by  the 
sli^t  reduc±ion  in- plasma  size  at  hi^ier  flew  rates. 

Die  average  temperatures  in  Fig.  33  eu?e  seen  to  be  cxinsiderably 
less  than  the  peak  temperatures.  Di^  range  from  about  2200  K  at  lew 
flow  rates  to  about  1200  K  at  hi^  flew  rates.  Further,  the  4  kW 
plasma  shews  a  sli^tly  hi^ier  mean  temperature  throu^iait.  Dius, 
better  performance  is  possible  at  lew  mass  flew  rates,  but  the  overall 
thermal  efficiencies  are  poor;  vhereas  at  hi^  mass  flew  rates, 
thental  efficiencies  are  hi^  but  the  performance  is  Icwer.  Clearly, 
a  cciiprcniise  has  to  be  made  between  these  two  opposing  requirements. 

A  similar  situaticai  exists  with  regard  to  power.  Hi^ier  pewers  mean 
better  performance  but,  as  seen  in  Figs,  31  and  32,  hi^  pewers  are 
also  less  efficient. 

Finally,  to  get  an  idea  of  vhat  level  of  performance  th^e 
temperatures  premise,  we  calculate  the  Igp  (specific  impulse)  using 
simple  one-dimensional  expansion  to  vacuum.  Corresponding  to  a 
temperature  of  2200  K,  the  I^  is  480  s  while  at  the  lower  temperature 
of  1200  K,  the  Igp  is  350  s.  Dhese  performaiK«  numbers  are  in  the 
range  that  current  chemical  systems  can  provide.  Clearly,  the 


feasibility  of  wan^guide  h^ted  prcpulsicn  ^  of  oacrciwave 
prqp^idn,  in  general,  depends  on  \iAiether  hig^hf^  p^onaanoe  can  be 
attzuned.  Diis  will  be  an  iii|xirtant  a^sect  of  inicxow^  pEcopulsion 
research  in  ^  next  fw  yee^. 
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4.0 

A  detailed  analytical  study  of  mcxowave  prcpulsicn  has  been 
undertaken.  Ihe  absorpticai  and  flcwfield  characteristics  of  three 
different  aicrcwave  absorption  coifigurations  —  free-floating  plasma 
in  resonant  cavities,  bluff-bocty  stabilized  plasmas  in  resonant 
cavities,  and  bluff-^xx^  stabilized  plasmas  in  waveguides  —  have  been 
inclixSed.  Results  have  been  ocnpared  with  experimental  measurements 
both  qualitatively  and  quantitatively'.  In  general,  the  model 
correctly  predicts  the  she^ie,  size  and  locatiai  of  the  plasna,  and  its 
response  to  increases  in  microwave  pcwer,  changes  in  the  pressure,  the 
plasma  movement  with  changes  in  the  flew  velocity,  and  the  expansion 
of  the  single  plasma  to  a  double  plasma  in  several  limiting 
ocxKditicns. 

In  terms  of  quantitative  ccnpariscais,  the  model  prc^jerly  predicts 
the  peak  tenperature  in  the  plasma  for  the  various  geometric 
configurations  tested.  Of  particular  interest  here  is  the  fact  that 
even  thou^  the  peak  temperature  is  nearly  constant  for  any  cxie 
ccsifiguraticxi,  the  model  correctly  predicts  the  differences  in  peak 
temperatures  that  have  been  observed  experimentally  between 
geometrical  configurations.  In  particular,  experimental  measurements 
of  the  peak  tenperatures  in  free-floating  plasmas  were  in  good 
agreement  with  previous  microwave  plasma  tenperature  measurements  by 
other  investigators,  but  the  measurements  in  bluff-bocty  stabilized 
plasmas  had  shewn  peak  tenperatures  that  were  considerably  hi^er  than 
these  well-established  levels.  Consequently,  it  had  been  thou^t  that 
th^  were  in  error. '  Uie  predictions,  however,  shew  that  the  peak 
tenperature  in  bluff-botty  stabilize  plasmas  increases  substantially 
as  the  bo<^  size  is  decreased,  and  predict  values  in  close  agreement 
with  the  experiments  for  the  same  sized  bluff  bodies.  In  addition  to 
the  taiperature  oenparisens,  the  predicted  coupling  efficiencies  are 
in  good  agreement  with  the  experimental  measurements  also.  Finally, 
predictions  of  the  low  power  threshold  for  naintcuning  a  plasma  are 
well-predicted  for  all  three  geonetries. 

Bie  predictions  indicate  that  the  experimental  conditions  have 
all  been  run  at  conditions  that  are  near  the  boundary  in  vAiich 
buoyancy-induced  recirculation  is  observed  in  the  plasma.  In  the  flow 
ranges  of  the  experiments,  snail  changes  can  lead  to  strong 
recirculatiai,  or  can  remove  the  recirculating  region  entirely, 
di^)ending  cxi  the  pressure  and  velociiy  levels.  Cenputations  with  and 
without  gravitatioal  effects  included  show,  hewever,  that  the  plasmas 
are  not  strongly  affected  by  this  recirculation,  so  that  the 
experimental  results  are  not  suspect  because  th^  are  conducted  near 
this  boundary.  The  analytical  results  also  show  that  the  presence  of 
the  plasma  with  its  strong  heating  characteristics  can  reduce  or  even 
remove  the  familiar  recirculation  region  that  occurs  behind  bluff 
bodies,  but  that  they  centinue  to  serve  effectively  as  "flame- 
holders". 
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A  major  shortxxming  of  the  micaxiwave  plasma  experiments  that  have 
been  mcide  to  date  is  that  althou^  they  have  all  shown  hi^  peak 
temperatures,  the  average  gas  temperatures  have  been  quite  low.  Ihese 
small  mean  temperature  rises,  seldcmi  more  than  a  few  hundred  Kelvins, 
strongly  limit  the  expected  qiecific  impulse  provided  by  potential 
microwave  thermal  thrusters.  Ihe  parametric  studies  based  on  the 
model  suggest  ways  for  increasing  average  temperature  to  values  that 
will  produce  attractive  ^jecific  impulses.  For  the  bluff-bocfy 
stabilized  plasmas,  the  average  temperature  increases  rapidly  as  the 
size  of  the  blxoff  bcx^  is  increased.  Ihe  reascxi  is  that  the  bluff- 
bocfy  size  sets  the  diameter  of  the  plasma,  as  well  as  the  amount  of 
mass  flow  that  goes  throuf^  the  plasma.  Similarly,  in  the  free- 
floating  plasma,  the  fraction  of  fluid  that  flows  throu^  the  plasma 
can  be  increased  by  carefully  CcntroUing  the  tube  diameter,  its 
shepe,  and  the  mass  flow.  Such  cptimiizations  are  best  made 
analytically  to  be  able  to  identic  the  regimes  in  vhich  ejqarimental 
demonstraticsns  shculd  be  made,  vhile  still  keying  the  plasma 
sufficiently  far  from  the  wall  to  prevent  arcing.  Additicmal  research 
in  this  area  is  needed  to  (demonstrate  noninal  ipper  bound  temperature 
measurements  that  can  be  obtained  with  micrcwave  plasmas  to  properly 
(define  their  practical  ranges  of  specific  impulse. 

The  otha:  major  limdtation  of  micrcwave  thrusters  is  the  maximum 
thrust  size  limiitation  set  by  the  wavelength  of  the  micrcwaves.  The 
wavelength  inherently  sets  the  size  of  the  plasma,  and,  hence,  the 
ncminal  thrust  level  because  absorption  is  most  effectively  limited  to 
resonant  cavity  absorpticxi.  Waveguicde  absorption  has  been 
demonstrated  herein,  but  it  is  not  as  efficient  (in  terms  of  cxxpling 
efficiencies)  as  rescxant  cavity  stabilized  plasmas.  Waveguicde 
plasmas  reflect  substantial  fractions  of  the  incident  pcwer  back  to 
the  load  vhere  they  are  lost  and,  in  addition,  they  premise  little  in 
terms  of  increased  thrust  levels.  Because  of  these  size  limitations, 
it  is  impossible  to  heat  large  quantities  of  gas  in  a  micrcwave 
thruster,  so  that  even  theu^  hi^  pewer  micrcwave  scurces  can  be 
built  in  li^twei^t  form,  they  cannot  be  used  (directly  to  scale  ip 
the  thrust  sizes  available  with  micrcwave  thrusters.  Specific  efforts 
directed  toward  this  scaleip)  process  have  never  been  attempted  and 
ocnc^Ttual  studies,  using  the  available  analytical  and  ejqperimental 
tools  edesexibed  herein,  would  appear  to  be  a  second  premising  area  for 
further  research. 
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Figure  la.  Fluid  dynamic  grid  for  sphere-cylindrical  tube  with  the 

microwave  cavity.  Right  side  of  the  cavity  shows  the  no-loss  field 
pattern. of  the  TM012  mode. 
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Figure  lb.  Geometry  for  bluff-body  stabilized 
resonant  cavity. 
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Figure  2.  Schematic  of  experimental  set-up  for  plasmas  generated  in  a  circular  waveguide 


Figure  3. Undistorted  axial  electric  field  in  a  waveguide,  (a)  waveguide  with  reflecting 
end  wall,  (b)  transmitting  waveguide  with  one- dimensional  method  of  charac¬ 
teristics,.  and  (c)  transmitting  waveguide  with  proper  characteristic  boundary 
condition. 


Figure  4.  Fluid  dynamic  grid  for  axisymmetric  bluff  body  geometry. 
Flow  is  from  top  to  bottom. 
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Figure  5 .  Convergence  of  low  Mach  nuxhber  computations  for  the  bluff 
body  "geometry  at  different  Reynolds  numbers. 


Figure  6.  Bluff  body  computations  for  different  Reynolds  numbers. 
Re  =  5,  60,  350, 1000. 
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Figure  7 .  Convergence  of  coupled  Navier-Stokes  and  Maxwell  equations  for  the 
bluff  body  geometry  using  the  TMqu  mode. 
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Figure  9 .  Effect  of  input  power  on  bluff  body  plasmas. 
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Figure  11.  Representative  solution  for  bluff-body  stabilized  ■ 
resonant  cavity  plasma. 


Figure  12.  Coupling  efficiency  of  bluff-body  stabilized 
cavity  plasmas.  Comparison  of  predictions 
and  experiments. 


re  14.  Representative  solution  of !the  sphere-cylinder  configuration. 
P  =  1  atm,  771  =  1  X  10"®  kg/s,  Re=10,  Pi„c  =  3  kW. 


o 


CNJ 

d 


il 


V 

2  O 

O 

Mtn 

:^’ 

^  V 

<i>  •>• 

^ 'd 

to  a) 

•e  S 
0-5 
'43^ 

s  s 

o  ^ 


to 

'O 

V 

cc 


o 

d 

0^ 

t-i 


L)  W- 

•d,  V 

U»  rj 

"5  :d 

u  .-^ 
_ 

T)  cr5 

«  . 

O 

(H  •-* 

Q 


.2  « 

Q '« 


LO 


u 

(-1 

3 

bO 


52 


STREAMLINES 

I  =;  0.02.  m/s  0-\l  m/s  O  '  5  m/s  1*1  m/s 

=  -2.  12.  50  Izo. 

Figure  19.  Effect  of  gas  velocity  variation  on  straight  duct  plasmas. 

Po  =  1  atm,  P  =  Z  icW. 


Figure  21.  Effect  of  discharge  pressure  on  straight  duct  plasmas. 
P  =  3  kW. 
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Figure  23.  Coupling  efficiency  of  free-floating 

cavity  plasmas.  Comparison  of  predictions 
and  experiments. 
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Figure  31.  Coupling-efficiency  as  a  function  of  mass  flow  rate. 


